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Non-compact QED-j, is simulated both in the quenched and unquenched cases. In particular, we investigate 
the restoration of chiral symmetry at finite temperature. We also compute the zero temperature spectrum of 
the theory, including (in the quenched case) the dynamical fermion mass. From these two set of data, one can 
obtain estimates for the ratio of the mass gap to the critical temperature, of particular interest for applications 
to high— Tc superconductivity. 



1. Introduction 

Over the last few years, QED^ has attracted a 
lot of attention because of potential applications 
to models of high— Tc superconductivity. How- 
ever, at the present time, the phase diagram of 
the theory is still not known precisely. There 
have been several analytical studies but approx- 
imations are often drastic and different results 
are obtained depending on what assumptions are 
made. Hence the interest of trying to see if a 
lattice gauge theory simulation can sort out this 
situation. 

To avoid any confusion, it is worth mention- 
ing that we are considering here the 4-component 
formulation of QED^ where the Dirac algebra is 
represented by the 4x4 matrices 70, 71 and 72. 
This formulation preserves parity and possesses 
an SU(2) group of chiral symmetries generated 
by 73, 75 and 7375. It also arises naturally as 
the continuum limit of the staggered fermion con- 
struction in (2+1) dimensions ||l|. If a mass term 
is spontaneously generated, the above SU(2) sym- 
metry will break down to the U(l) generated by 
7375 and there will appear two Goldstone bosons 
associated with the operators V'73'0 and ip"/5ip. 

Before discusssing the numerical simulation, it 
is useful to make a brief review of the existing 
analytical results. These are all based on approx- 
imate treatments of the Schwinger-Dyson equa- 
tion for the fermion self-energy. The first approx- 
imation which was considered is the ^/ Nf expan- 
sion {e^Nf fixed) where Nf is the number of 



fermion flavors. To first order, there is no vertex 
correction and the photon propagator becomes: 
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A numerical solution of the Schwinger-Dyson 
equation then gives for the critical number of fla- 
vors beyond which chiral symmetry is restored: 
Nc = 3.28. This result has been criticized by Pen- 
nington and Walsh |^] , who by considering a plau- 
sible expression for the interaction vertex arrive 
at the conclusion that chiral symmetry is always 
broken (however high Nf can be). The finite tem- 
perature theory has also been investigated Q . In 
this case, the longitudinal photons are screened, 
but the transverse photons remain unscreened: 
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The transverse photons therefore give rise to se- 
vere infrared divergences and it is not known at 
present how this problem can be cured. Dorey 
and Mavromatos have assumed that transverse 
photons can be neglected in the Schwinger-Dyson 
equation and have computed the critical temper- 
ature in this case. They obtain in this way for 
Nf = 2 a ratio of the mass gap to the critical tem- 
perature which is of order 10 (much bigger than in 
BCS-like theories, hence the potential interest for 
models of high- Tc superconductivity). Because of 
all the uncertainties in the above mentioned re- 
sults, it would be very interesting to have lattice 
measurements of the relevant quantities. We will 
first discuss the techniques that we have at our 



2 



disposal to carry out such a project and will then 
present our current results for the chiral conden- 
sate (both at zero and non-zero temperature) and 
for the meson spectrum. 

2. Symmetry Breaking on the lattice 

When studying chiral symmetry on the lat- 
tice, we encounter from the start a fundamental 
problem: since the number of degrees of freedom 
is finite, a symmetry can not be spontaneously 
broken. Therefore (iptp) will tend to zero when 
m —^ 0. However, one can hope that there will 
be two regimes in the dependence of {'ip'4') on m. 
In the intermediate mass range, (ipip) would be 
essentially independent of the size of the lattice, 
whereas for small masses, finite-size effects would 
show up and (V'V') would drop sharply to zero. 
If this is the case, then we can obtain an es- 
timate of the infinite volume value of (ip'^P) by 
extrapolating from the intermediate mass range. 
Although this technique usually provides useful 
information, it may not be powerful enough to 
really identify a phase transition. However, fur- 
ther checks are available to confirm the presence 
or the absence of a phase transition. To be sure 
that chiral symmetry is really broken in the con- 
tinuum, one has to show that {^ip) obeys the re- 
quired scaling behaviour with respect to (3. In 3 
dimensions (3^{4!ijj) (= continuum value of (V'V') 
in units of e^) has to be constant. In practice, one 
is looking for a range of f3 where there is a plateau 
in /3'^{4>ip). Note however that finite size effects 
can significantly delay the appearance of such a 
plateau. In Q for example it was shown that chi- 
ral symmetry is indeed broken in quenched QED3 
but lattices of size ranging up to (80)'^ were nec- 
essary for that purpose. In the case where the 
extrapolation method seems to indicate the pres- 
ence of a phase transition at finite /3, we can as- 
sess its existence by using techniques borrowed 
from the "theory" of critical phenomena. At the 
critical value of the coupling, for example, it is 
expected that {^ip) will behave like m^^^ (where 
d is one of the critical exponents characterizing 
the transition). 

Independently of these techniques based on 
measurements of the chiral condensate, one can 



also gain information about the phase diagram 
by considering the meson spectrum. One has the 
following alternative: Either there are parity dou- 
blets in the limit of massless fermions (and the tt 
and the a would be degenerate for example) or 
this is not the case and the mass of the pion sat- 
isfies the usual PCAC relation. Consequently for 
m=0, the ratio M^/Mq- should be a step function 
centered at /3 = /3c- If m is not exactly zero then 
the function will be smoothed. This has recently 
been given a quantitative form by applying the 
hypothesis of correlation length scaling . 

3. The chiral condensate at T=0 

We have extended the results of Q by measur- 
ing the chiral condensate on a (12)3 lattice. Put 
together, these results indicate that finite size ef- 
fects are strong, which makes the application of 
the techniques mentioned above rather difficult. 
But, if we try anyway, the data presented in fig.l 
amd 2 seem to indicate the presence of a phase 
transition around (3 = 0.2 for Nf = 4 and around 
f3 = 0.26 for Nf = 3 (although the later is less 
compelling because several curves are compatible 
with linear behaviour). The critical exponents 
would be respectively 6 = 3.7 and 2.7 but the un- 
certainties are large. Clearly these results require 
some further confirmation either from measure- 
ments on bigger lattices or from some indepen- 
dent determination of S. 

4. Spectrum computations 

We used Kogut-Susskind fermions and com- 
puted the correlators between mesonic operators 
which are local in the lattice variables. Because 
of lack of space, only the results for the pion will 
be presented here. We have accumulated data 
for Nf = 0, 2 and 4 at values of the coupling 
ranging from (3 = 0.2 to (3 — 0.4 and for masses 
m=0.0125, 0.025 and 0.05. In the quenched case, 
chiral symmetry is clearly broken: we observe the 
behaviour expected from PCAC over a wide range 
of values of (3 (Fig. 3) and the finite-size effects 
can be investigated systematically. For Nf = 2 
(Fig.4), at /? = 0.2, we observe on a (12)^ x 24 lat- 
tice that chiral symmetry is broken. At (3 = 0.4, 
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strong finite-size effects are present but a com- 
parison of our (12)^ X 24 and (14)^ x 24 results 
suggest tfiat cfiiral symmetry will still be broken 
in the infinite volume limit. For Nf = 4 (Fig. 5), 
chiral symmetry is broken at (3 = 0.2 but prob- 
ably not at (3 = 0.4 although finite size effects 
remain important even on a (16)^ x 24 lattice. If 
the above results are correct, they would confirm 
the analytical expectations based on a 1/iV/ ex- 
pansions, namely that the critical value of Nf be- 
yond which chiral symmetry is restored is around 
3. At this point, it should be mentioned that 
in their lattice investigation of compact and non- 
compact QED3, Azcoiti and Luo |^ arrive at a 
rather different picture of the phase diagram than 
the one suggested here. 

5. The chiral condensate at T ^ 

Our measurements of {i^ijj) at Nf — 2 on 
(12)2 ^ 4 j^jj^ (-g)2 X 4 lattices are shown on Fig. 6. 
Using the extrapolation method mentioned pre- 
viously, we would conclude that chiral symme- 
try is broken at /? = 0.2 (Note that the extrap- 
olated value of (■0V') increases with the spatial 
size of the lattice) whereas it is not broken for 
/3 > 0.3. However, one has to be very careful 
at this point. Indeed, if this is really indica- 
tive of a transition in the continuum, then we 
should see again the transition on lattices with 
Nt = 8 and it should happen between /3 = 0.4 
and f3 — 0.6. But nothing like this is observed on 
Fig. 7, instead all the curves extrapolate smoothly 
to zero. In order to further test for the presence 
of a transition, we can plot In(V'V') versus In to 
(Fig. 8) and see if there is a curve which exhibits 
critical behaviour. Only /3 = 0.3 would be com- 
patible with a straight line provided that we elim- 
inate the point at to = 0.00625. The presence of 
this point in fact seems to indicate that P = 0.3 
is an upper bound for the value of the critical 
temperature. However, one should consider the 
possibility that this point is sensitive to finite- 
size effects (although for (3 > 0.4 further runs on 
(24)2 ^ g lattices indicated that finite-size effects 
are small) and it is tempting to speculate that 
the actual transition is not far from /3 = 0.3. If 
this were the case, the ratio of the critical temper- 



ature to the square of the electric charge would 
be close to T/e^ = 0.0375, almost an order of 
magnitude bigger than the result of Dorey and 
Mavromatos (0.0043)! Further runs at lower val- 
ues of f3 are needed in order to resolve the am- 
biguity. Ideally one would like to obtain a curve 
with positive curvature on Fig. 8 since it would 
give us a lower bound on the value of the critical 
temperature. Finally, measurements of screening 
length (analogous to meson spectrum calculations 
at zero temperature) could also provide useful in- 
formation for locating the position of the phase 
transition. We hope to be able to report on this 
in the near future. 

6. Conclusion 

Our results seem to give some support to the 
following two statements: 

1. At zero temperature, chiral symmetry is re- 
stored beyond some critical value (iVc) of 
the number of flavors and Nc is close to 3. 

2. For Nf = 2, the temperature at which chi- 
ral symmetry is restored is much higher 
than the current analytical prediction Q. 

However, as is clear from the discussion given in 
the text, further work is still needed in order to 
completely clarify the situation. 
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